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摘要 :【 目 的】 分 析 中 国 冰 清 绢 蝶 Parnassius glacialis 遗传 多 样 性 和 遗传 分 化 情况 及 其 系统 发 生 关 
A ,推测 其 起 源 及 分 此 时 间 , 并 探讨 其 历史 生物 地 理 分 布 格局 的 成 因 。【 方 法 】 以 2008 -2016 年 采 
集 于 中 国 13 个 地 区 的 冰 清 绢 蝶 的 325 头 及 绢 蝶 属 其 他 11 个 种 的 11 头 个 体 样 品 , 对 其 线粒体 基因 
组 的 AT 富 集 区 (控制 区 ) 序 列 进行 PCR 扩 增 和 序列 测定 ,采用 MEGA 6.0, Dna SP 5.1 和 Arlequin 
3.5 等 软件 分 析 其 遗传 多 样 性 和 遗传 分 化 情况 ; 以 其 他 近 缘 绢 蝶 种 类 作为 外 类 群 ,采用 PhyMI3. 0, 
MrBayes 3.2 和 BEAST V1.8.3 等 软件 重建 冰 清 绢 蝶 的 系统 发 生 树 , 并 利用 宽松 分 子 钟 和 校准 方法 
推测 中 国 冰 清 绢 蝶 的 起 源 和 分 歧 时 间 ; 同 时 结合 现今 冰 清 绢 蝶 的 地 理 分 布 特 点 和 第 四 纪 以 来 的 地 
球 环境 事件 ,分 析 其 扩散 路 径 ,探讨 其 历史 生物 地 理 分 布 格局 及 成 因 。 【结果 】 冰 清 绢 蝶 AT 富 集 区 
序列 长 度 介 于 487 ~495 bp 之 间 , 平 均 为 491 bp ,其 长 度 的 差异 主要 体现 在 Poly(T) 或 Poly(A) 的 
数目 差异 ,A+T 平 均 含 量 为 935.76% 。 供 试 冰 清 绢 蝶 13 个 地 理 居 群 325 头 个体 中 , 共 检 测 出 基于 
线粒体 基因 组 AT 富 集 区 序列 的 239 个 单 倍 型 , 单 倍 型 多 样 性 Hd 为 0.9971; X ER 2 ASTE d8 a 
值 为 0.02948 ，Theta (per site) Eta 值 为 0.06594。 系 统 发 生 分 析 和 分 子 定年 显示 , 冰 清 绢 蝶 和 和 白 
绢 蝶 P. stubbendorfü 为 近 缘 姊妹 种 ,二 者 在 距 今 7.52 百 万 年 ( Ma) 的 中 新 世上 晚期 开始 分 化 。 冰 清 
绢 蝶 的 祖先 可 能 起 源 于 我 国 西南 地 区 (现今 横断 山脉 -喜马拉雅 附近 的 青藏 高 原 东北 缘 ) 距 今 1.53 
Ma 的 更 新 世 时 期 , 先 由 起 源 地 先期 扩散 到 达 小 陇 山 秦岭 一 带 ; 后 来 ,由 于 第 四 纪 的 冰期 - 间 冰 期 轮 
回 事件 和 扩散 路 径 的 不 同 而 分 化 出 2 大 支 系 , 但 两 个 支 系 之 间 存 在 少量 的 单 倍 型 相互 混杂 现象 ; 随 
后 ,2 大 支 系 继续 向 低 海拔 、 低 纬度 地 区 扩散 , 即 通过 伏 牛 山脉 向 东北 方向 扩散 到 党 山 泰山、 昆 党 
山 一 带 , 经 大 别 山系 向 东南 方向 扩散 到 了 琅 划 山 、 紫 金山 和 天 目 山 一 带 。【 结 论 ] 中 国 13 个 居 群 的 冰 
清 绢 蝶 没 有 形成 明显 的 地 理 分 化 格局 ,其 遗传 变异 主要 来 自居 群 内 部 , 居 群 间 的 遗传 分 化 现象 尚 不 
明显 ; 冰 清 绢 蝶 起 源 于 更 新 世 时 期 的 较 高 海拔 地 区 ,由 于 第 四 纪 的 冰期 - 间 冰 期 轮回 事件 以 及 扩散 
路 径 的 不 同 而 形成 2 大 支 系 并 向 中 国 东部 和 南部 低 海拔 和 低 纬度 地 区 扩散 。 

关键 词 : ERA. 冰 清 绢 蝶 ; 线粒体 基因 组 ; AT 富 集 区 ; 谱系 发 育 ; 遗传 分 化 ; 谱系 生物 地 理 
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Abstract: [ Aim] This study aims to analyze the genetic diversity, genetic differentiation, and 
phylogenetic relationships of populations of the alpine butterfly Parnassius glacialis in China, to infer their 
origin and divergence time, and to explore their spatiotemporally evolutionary history. [ Methods] The 
mtDNA AT-rich sequences of 325 individuals from 13 geographic populations of P. glacialis in China and 
11 individuals of other 11 Parnassius species were amplified and sequenced, and their sequence variation 
and genetic differentiation were analyzed using MEGA 6.0, Dna SP 5.1, Arlequin 3. 5 and other genetics 
analysis software. Using the closest relatives as the outgroups, the P. glacialis phylogenetic trees were 
reconstructed with such analytical software as PhyML3. 0, MrBayes 3. 2, and BEAST V1.8.3, and the 
origin and divergence time of the P. glacialis lineages were estimated by using relaxed molecular dating 
method and calibrations of the previous studies. Based on the present biogeographic distribution of P. 
glacialis and the main earth. environmental events since the Quaternary Period, its dispersal route and 
spatio-temporal pattern. and the underlying earth environmental factors were tentatively inferred. 
[ Results] The AT-rich sequences of 325 P. glacialis individuals from 13 geographic populations ranged 
from 487 to 495 bp in length with the mean of 491 bp, being relatively conserved in length. The main 
length difference of the AT-rich region lies in the Poly-T or Ploy-A copy number, and the mean A +T 
content of the region is up to 95. 7696. Totally 239 halotypes based on mtDNA AT-rich sequences were 
found in the 325 individuals of the 13 populations of P. glacialis, with its haplotype diversity ( Hd) value 
of 0. 9971, the nucleotide polymorphism index « of 0. 02948, and the Theta (per site) Eta value of 
0. 06594. The reconstructed phylogenetic trees and molecular dating showed that P. glacialis is most 
closely related to P. stubbendorfü, and these two sister species began to differentiate at about 7. 49 
million years ( Ma) ago in the late Miocene. The ancestors of P. glacialis were originated in the 
southwest of China (near the Hengduan Mountain-Himalaya Mountain of the Qinghai-Tibet northeastern 
edge nowadays) at about 1. 53 Ma ago; then, they dispersed from their originating site into Xiaolongshan 
and Qinling Mountains; later, they diverged into two main lineages with some inter-mingled halotypes 
along the mountainous routes, accompanied by the Quaternary glacial-interglacial changes; lastly, the 
one lineage continuously spread northeastward to Funiu, Songshan, Taishan, Kunyushan Mountains, and 
the other lineage spread step-by-step to Dabie, Langyashan, Zijinshan and Tianmushan Mountains. 
[Conclusion] The results of our phylogenetic analysis of P. glacialis indicate that no pronounced 
correlations exist between genetic differentiation and geographical distance judged from higher genetic 
variance within populations than between populations. In addition, we propose that the divergence of the 
two main lineages of P. glacialis occurred in the mountainous areas with relatively higher latitude in the 
Middle Pleistocene and then spread to the mountainous areas with relatively lower latitude in South and 
East China due to the Quaternary glacial-interglacial cycles and different dispersal routes. 
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2H UE iS Parnassius 隶属 于 风暴 科 ( Papilionidae ) 
绢 蝶 亚 科 ( Parnassinae ) , © MTER rp — A F2) 38 
特 的 类 群 ,其 中 的 绝 大 多 数 种 类 生活 在 高 海拔 的 高 
原 地 带 ,特别 是 我 国 的 青藏 高 原 地 区 。 由 于 长 期 对 
特殊 地 质 .气候 等 环境 的 适应 性 进化 ,它们 在 形态 和 
生理 上 逐渐 产生 了 一 系列 高 海拔 环境 的 适应 性 特 
征 ,如 耐 低 温 、 耐 缺 氧 .和 抗 紫外 线 等 ( 周 范 ，1999 ; 





























武 春生 和 徐 博 峰 , 2017) 。 目 前 ,它们 已 成 为 昆虫 物 
种 形成 与 演化 .适应 性 分 子 进化 等 研究 领域 十 分 难 
得 的 模式 生物 类 群 之 一 ( 夏 琛 琛 , 2016; Carvalho et 
al., 2017; Rytteri et al., 2017; E3 R, 2017; 张 敏 
等 , 2017) 。 前 期 的 有 关 研究 表 明 , 绢 蝶 属 祖先 很 可 
能 起 源 于 第 四 纪 早 期 ( 距 今 20 Ma 左右 ) 的 我 国 西 
南部 地 区 (现今 青藏 高 原 周边 一 带 ) ( 刘 文 萍 ， 
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1999; Tzedakis et al., 2002; Michel et al., 2008; £ 
运 良 , 2017) ,之 后 , 随 着 地 球 环境 的 不 断 演变 和 自 
身 的 快速 适应 辐射 而 向 其 他 地 区 扩散 。 其 中 , 冰 清 
绢 蝶 Parnassius glacialis 是 我 国境 内 分 布 海拔 相对 
较 低 、 分 布地 域 相对 广泛 ( 西北、 华中 和 华东 的 绝 大 
部 分 地 区 ) 、 也 是 唯一 分 布 到 我 国 江 南 地 区 的 一 种 
ARBER (Ji] 3e , 1999 ; 武 春生 和 人 徐 博 峰 , 2017) 。 
迄今 ,关于 其 起 源 分 化 的 时 空格 局 及 相关 的 地 球 环 
境 背 景 尚 不 清楚 ,值得 深入 探讨 。 

昆虫 的 线粒体 基因 组 通常 是 一 类 15 ~16 kb 的 
双 链 环 状 DNA 分 子 , 具 有 分 子 量 小 、 结 构 简 单 .组 成 
稳定 .母系 遗传 和 进化 速度 快 等 特点 , 目前 已 成 为 探 
讨 生物 类 群 较 低 等 级 分 类 阶 元 间 系 统 演化 的 理想 工 
具 之 一 ( Wolstenholme, 1992; Boore, 1999; Nei and 
Kumar, 2000; Simon et al., 2006) 。 其 中 ,AT 富 集 
区 (AT-rich region ) 或 控制 区 (control region ) 是 整个 
线粒体 基因 组 序列 和 长 度 变 异 最 大 的 区 域 , 属 非 编 
码 区 ,主要 参与 线粒体 DNA 复制 和 转录 的 调控 。 它 
们 在 进化 过 程 中 所 受 的 选择 压力 较 小 ,进化 速率 快 ， 
含有 较为 丰富 的 多 态 性 。 目 前 已 被 广泛 应 用 于 物种 
内 部 的 分 子 系统 学 、 分 子 进化 学 .种群 遗 传 学 和 谐 系 
生物 地 理学 等 研究 领域 ( 郭 松 长 , 2006; 朱文 奇 等 ， 
2010; JK K BUR TL tE $, 2011; Huang and Ke, 
2014; Li et al., 2014; 要 晓 旭 等 , 2015; 何 海燕 等 ， 
2016; 涂 世 英 等 , 2016; 章 明 等 , 2016 ) 。 

本 人 研究 尝试 对 采集 于 我 国 西北 、 华 中 以 及 华东 
等 不 同 地 域 的 13 个 冰 清 绢 蝶 地 理 居 群 共 325 头 个 
体 的 线粒体 基因 组 AT 富 集 区 序列 进行 PCR 扩 增 、 
序列 测定 和 比较 分 析 ; 基 于 获得 的 序列 数据 ,运用 多 
种 遗传 分 析 方法 估算 群体 间 和 群体 内 的 遗传 分 化 状 
况 ; 同 时 ,基于 不 同 序列 进化 模型 和 多 种 建树 方法 重 
建 它们 的 系统 发 生 树 ,依据 前 期 的 分 子 定年 研究 结 
果 以 宽松 分 子 钟 方法 估算 各 居 群 间 的 分 于 时 间 ,并 
结合 第 四 纪 以 来 我 国 的 主要 地 球 环境 变化 事件 等 地 
史 资 料 ,初步 探讨 冰 清 绢 屿 在 我 国 的 起 源 .谱系 生物 
地 理 分 化 格局 及 其 地 球 环境 成 因 。 

























































































1 材料 与 方法 


1.1 实验 材料 
冰 清 绢 蝶 采 集 时 间 为 2008 -2016 年 ,采集 地 见 





1 -25 ,用 于 实验 分 析 。 作 为 外 群 同样 用 于 实验 分 析 
的 绢 蝶 属 其 他 11 个 种 的 11 头 个 体 标 本 的 采集 信息 
未 列 出 。 
1.2 冰 清 绢 蝶 线 粒 体 基 因 组 AT 富 集 区 序列 扩 增 
和 测定 

使 用 动物 基因 组 快速 抽 提 试剂 盒 ( 上 海 生 工 生 
物 工程 有 限 公 司 , 上 海 ) ,参照 使 用 手册 提取 冰 清 绢 
蝶 总 基因 组 DNA。 参 考 Simon 等 (2006 ) 并 结合 冰 
清 绢 蝶 12S rRNA 和 ND2 的 基因 序列 保守 区 进行 引 
物 设 计 ( 王 运 良 , 2017 ;线粒体 全 基因 组 序列 已 测 
出 ,尚未 提交 至 GenBank): CRIF: 5'-ATAATAGGG 
TATCTAATCCTAGT-3'; CRIR: 5'-CCAATTCAACAT 
CCCAATCA-3', CR2F; 5'-ATAGGGTATCTAATCCT 
AG-3'; CR2R; 5'-CTAATCCAATTCAACATCC-3', 
用 于 扩 增 128 rRNA , AT 富 集 区 tRNA“ 和 tRNA™ 
的 一 部 分 ,产物 长 度 约 为 1.2 kb。 引 物 由 上 海 生 工 
生物 工程 有 限 公 司 合 成 。PCR 反应 体系 为 50 uL: 
Taq DNA 聚合 酶 (5 U/uL)0. 6 uL, Buffer (KE 
Mg )6.0 uL, MgCl, (25 mmol/L) 8. 0 uL, dNTPs 
(2.5 mmol/L)1.5 pL， 上 下 游 引物 (10 pmol/L) 4& 
1.8 uL, DNA 模板 (1.0 ng/uL) l. 5 uL, ddH,0 
28.8 uL, PCR 反应 条 件 : 95Y 预 变性 5 min; 然后 
95% 变性 40 s, 52% 退火 1 min, 72% 延伸 2 min , 共 
35 个 循环 ;再 行 72%C 总 延伸 10 min。PCR 反应 所 用 
试剂 均 购 于 上 海 生 工 生物 工程 有 限 公 司 ,PCR 扩 增 
产物 经 分 离 纯化 后 由 华 大 基因 合肥 世 诚 生物 有 限 公 
司 进 行 双 向 测序 。 
1.3 序列 处 理 

将 1.2 节 测 序 所 得 结果 仔细 与 测序 峰 图 核对 校 
正 ,使 用 Clustal X 2.1 (Larkin et al., 2007 ) 对 测序 
结果 进行 比 对 后 拼接 ,然后 在 NCBI. 上 使 用 Blast 进 
行 序列 同 源 性 搜索 ,以 验证 所 得 序列 是 否 为 目标 片 
段 。 拼 接 后 的 序列 以 Clustal X 2.1 进行 序列 多 重 比 
对 和 辅 以 人 工 校 正 ,应 用 MEGA 6.0 (Tamura et al., 
2013) 进行 碱 基 舍 量 以 及 多 态 位 点 分 析 , 并 基于 
Kimura 双 参 数 ( Kimura 2-parameter, K-2P) 模型 计 
算 居 群 间 及 居 群 内 的 遗传 距离 ,利用 DnaSP 5. 10 
( Librado and Rozas, 2009) 进行 居 群 间 及 居 群 内 的 
遗传 多 样 性 分 析 。 
1.4 遗传 距离 及 多 样 性 分 析 

基于 DnaSP 5.10 统计 分 析 冰 清 绢 蝶 13 个 地 理 










































































Klo 标本 采集 后 立即 放 入 无 水 乙醇 中 国定 , 带 回 








居 群 325 头 个 体 的 基于 线粒体 AT 富 集 区 序列 单 倍 








实验 室 放 入 -80%C 冰 箱 保 存 。 本 研究 在 每 个 地 理 居 
群 所 采集 的 样本 中 随机 挑 取 各 25 头 个 体 ,分 别 编号 





型 个 数 、 单 倍 型 多 样 性 (haplotype diversity, Hd) 以 
及 核 背 酸 多 态 性 指数 普 等 遗传 多 样 性 指标 .应 用 
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表 1 样本 信息 及 基于 线粒体 基因 组 AT 富 集 区 序列 的 冰 清 绢 蝶 13 个 地 理 居 群 的 
单 倍 型 . 单 倍 型 多 样 性 和 核 背 酸 多 态 性 指数 
Table 1 Sampling information and the estimated haplotype, haplotype diversity and nucleotide diversity of 
13 geographic populations of Parnassius glacialis based on mitochondrial AT-rich sequences 
au 居 群 代码 ”样本 数 地 理 坐 标 a 单 倍 型 多 样 性 ，” 核 音 酸 多 态 性 指数 
样本 采集 地 . 单 倍 型 ; 
. . Population Number of Geographic Haplotype Nucleotide polymorphism 
Populations locality j Haplotypes 
code samples coordinates diversity Hd index ar 
YT Meus ye AM 
江苏 南京 紫金 山 ZJS 25 118. 83?E, 32. 06°N 19 0.953 0.02755 
Zijinshan, Nanjing, Jiangsu 
江苏 连云港 云 台 
江苏 连云港 云 合 山 YTS 25 119.40°E, 34.71°N 13 0.940 0.01906 
Yuntaishan, Lianyungang, Jiangsu 
nnd o2 gate 
安徽 六 安 天 党 于 TTZ 25 115. 77?E, 31. 17?N 20 0.965 0.03856 
Tiantangzhai, Liu'an, Anhui 
西 宝 鸡 黄柏 源 
i SRAI HBY 25 107. 40° E, 33.73°N 20 0.973 0.02777 
Huangbaiyuan, Baoji, Shaanxi 
河南 : 
NIBHÉFE AERE t LJS 25 111.66°E, 33.76°N 20 0.989 0.01505 
Laojunshan, Luoyang, Henan 
b Yr 
浙江 临安 天 目 山 TMS 25 119.45? E, 30.34°N 16 0.985 0.02746 
Tianmushan, Lin'an, Zhejiang 
东 泰 安泰 
Mice dun TS 25 117. 12?E, 36.25?N 16 0.983 0.02870 
Taishan, Tai'an, Shandong 
M 小 陇 
甘肃 天 水 小 陀 山 XLS 25 105. 68? E, 34. 85?N 23 0.985 0.03873 
Xiaolongshan, Tianshui, Gansu 
安徽 滁州 
SCENES LYS 25 118.29? E, 32.28?N 18 0.980 0.01844 
Langyashan, Chuzhou, Anhui 
MHRS 
HET 关山 NTS 25 110.75°E, 32.62°N 22 0.984 0.03066 
Niutoushan, Shiyan, Hubei 
ARA E Eb dg 
URAR ER KYS 25 121.73°E, 37.28°N 15 0. 983 0.04107 
Kunyushan, Yantai, Shandong 
ij 1 dn 
四 北 神农 染 SNJ 25 110.35?E, 31.52?N 21 0.967 0.04039 
Shennongjia, Hubei 
河南 郑州 涯 
河南 郑州 凯 山 SS 23 113. 05?E, 34.49°N 16 0.985 0.04199 


Songshan, Zhengzhou, Henan 


Arlequin 3. 5 ( Excoffier and Lischer, 2010) 中 分 子 变 
54r (analysis of molecular variance, AMOVA) 以 
估算 地 理 居 群 的 遗传 分 化 指数 (F-statistics，Fst) 以 
及 居 群 间 每 世代 的 迁移 数 (number of migrants per 
generation, Nm) ; 同时 ,基于 获得 的 控制 区 单 倍 型 序 
列 ,采用 DAMBE 7.0 (Xia, 2018) 把 其 序列 fas 格式 
转换 为 phy 格式 ,然后 通过 NetWork 4. 5. 1. (Posada 
and Crandall, 2001) 基于 中 间 连 接 法 (median 
joining) TJ ££ & sc [El , 
1.5 系统 发 育 分 析 

本 研究 将 获得 的 13 个 地 理 居 群 325 条 冰 清 绢 
WEB 12S rRNA 部 分 序列 、AT 富 集 区 tRNA” 和 
tRNA" 部 分 序列 进行 进行 序列 比 对 ,删除 AT 富 集 
区 中 不 稳定 的 高 变 区 ,得 到 冰 清 绢 蝶 各 居 群 的 单 倍 





















































型 序列 ,然后 以 白 绢 蝶 P. stubbendorfii 珍珠 绢 蝶 P. 
orleans 等 绢 蝶 属 其 他 11 个 种 作为 外 群 ,构建 单 倍 型 
序列 数据 集 Data 1。 运 用 最 大 似 然 法 ( maximum 
likelihood, ML) 和 贝 叶 斯 法 ( Bayesian inference, BI) 
构建 系统 发 生 树 。 最 大 似 然 法 分 析 时 , 先 使 用 软件 
PAUP* 4.0b10 ( Swofford, 2002 ) 和 ModelTest 
( Darriba et al., 2012) 联合 确定 其 最 适 进化 模型 
(GTR +I +G), Hii ii PhyML3.0 (Guindon et al., 
2010) 软 件 以 非 参 数 自 展 法 重复 抽样 1 000 次 ,获得 
系统 树 的 拓扑 结构 和 各 节点 的 自 展 检验 值 
( bootstrapping value) ; 贝 叶 斯 分 析 时 , 先 将 数据 集 的 
文件 格式 转变 为 nex 文件 格式 并 复制 到 软件 
MrBayes 3. 2 ( Huelsenbeck and Ronquist, 2001; 
Ronquist et al., 2012) 文件 夹 中 ,根据 Modeltest 估算 







































































NI 


4 HH Ej 














良 等 : 基于 线粒体 基因 组 控制 区 序列 的 





Pp 国 冰 清 绢 蝶 遗 传 分 化 和 谱系 生物 地 理 研究 





479 








的 最 优 碱 基 蔡 代 模 型 (GTR +1+ 6) 设置 参数 , 即 
lset 设置 替换 模型 nst =6( GTR 模型 ) ,位 点 速率 变 
异 模型 设置 为 rates = gamma ( Gamma 分 布 ) ,然后 运 
行 4 条 MCMC 链 ,以 随机 树 为 起 始 树 , 共 运 行 100 
万 代 , 每 100 代 抽 样 1 次 ,重复 1 次 ,运算 结果 显示 
分 列 频率 平均 标准 差 小 于 0. 01 时 ,舍弃 2 500 RE 
化 样本 ,剩余 样本 用 于 构建 系统 一 致 树 , 如 果 平 均 标 
准 差 没 有 小 于 0. 01 ,那么 增加 代数 ,继续 运行 。 根 
据 马 尔 科 夫 链 的 蒙特 卡 洛 方法 计算 各 节点 的 后 验 概 
率 值 (posterior probability ) 。 
1.6 谱系 发 育 年 代 学 分 析 

依据 系统 发 育 分 析 过 程 中 较为 稳健 的 BI 树 , 运 
JH BEAST V1.8.3 (Drummond and Rambaut, 2007 ) 
软件 ,基于 宽松 分 子 钟 模型 ,使 用 贝 叶 斯 算法 运算 构 
建 冰 清 绢 蝶 系 统 发 生 过 程 中 各 不 同居 群 分 支 事件 的 
时 间 表 。 在 分 析 过 程 中 ,结合 本 研究 前 期 的 研究 结 
果 ( 王 运 良 , 2017), 以 绢 蝶 属 的 分 化 时 间 (19. 81 
Ma, CI; 24.2 -15.4 Ma) 设 为 系统 树 根 部 的 时 间 界 
限 ,作为 校正 点 中 ,以 白 绢 蝶 和 珍珠 绢 蝶 的 分 歧 时 间 
(12.24 Ma) 为 校正 点 @, 并 以 白 绢 蝶 和 冰 清 绢 蝶 的 
分 下 时 间 (7. 49 Ma ) 作为 分 歧 时 间 推 测 的 校正 点 
,把 上 述 3 个 时 间 校 正点 嵌入 本 研究 BI 树 相应 的 
节点 位 置 ,运行 10 000 000 代 , 生 成 xml 文件 。 采 用 
Tracer 1.5 ( Rambaut et al., 2014) 打开 log. txt 文件 
评估 运行 结果 , 若 有 效 取 样 数 ESS ( effective sample 
size) (H 2:200 ,说 明 佑 算 结 果 比 较 理 想 , 否则 ,增加 
运行 代数 或 者 重新 调整 参数 继续 运行 。 运 行 结束 
后 ,将 单 次 或 多 次 生成 的 fas. txt 格式 文件 添加 到 


































































































BEAST v1.8.3 的 子 程序 Log Combiner vl. 8.3 rp yt 
行 合并 ,导出 合并 的 文件 (fas. tree 格式 ) ,运行 Tree 
annotator v1.8.3 把 之 前 导出 的 fas. tree 文件 转换 为 
mcc. trees 格式 文件 ;最 后 运用 FigTree vl. 4. 0 
(http: // tree. bio. ed. a-c. uk/software/figtree/) 打开 
BEAST 运行 后 生成 的 mec. tree 文件 ,获得 绢 蝶 属 各 
文系 及 冰 清 绢 蝶 不 同居 群 的 分 上 疏 时 间 。 




















2 结果 


2.1 冰 清 绢 蝶 线 粒 体 基因 组 AT 富 集 区 的 序列 组 
成 及 结构 特征 

本 研究 13 个 地 理 居 群 325 头 冰 清 绢 蝶 的 线 粒 
体 AT 富 集 区 序列 经 校对 后 提交 至 GenBank rp, X* 
录 号 为 MH561389 - MH561713。AT 富 集 区 序列 长 
度 介 于 487 ~495 bp 之 间 ( 平 均 为 491 bp) ,其 长 度 
差异 主要 体现 在 Poly (T). 或 Poly( A) 的 数目 差异 。 
T, C, A 和 G4 种 碱 基 的 平均 含量 分 别 为 47.7%， 
2.896, 48. 1% 和 1. 4%,A+T 平 均 含量 高 达 
95. 7696 ( 表 2)。 对 325 头 个 体 的 线粒体 DNA. AT 
富 集 区 序列 比较 分 析 发 现 , 所 有 冰 清 绢 蝶 的 AT 富 
集 区 序列 可 以 大 体 分 为 5' 端 保守 区 、 高 变 区 、3’' 端 保 
守 区 。5' 端 保守 区 在 tRNA 上游, 长 约 244 bp, HE 
中 包含 一 个 典型 的 约 11 bp 的 Poly( A) 结 构 , 它 们 在 
所 有 蝶 类 中 均 非 常 保守 。3' 端 保守 区 在 12S rRNA 
下 游 ,长 度 约 44 bp, 此 保守 区 包含 一 个 “ATAGA 
(T)” 7018 ~20 bp 的 Poly(T) 元 件 ,这 一 结构 在 绝 
大 多 数 蝶 类 中 也 广泛 存在 。 


























R2 13 个 不 同 地 理 居 群 的 冰 清 绢 蝶 线 粒 体 DNA AT 富 集 区 序列 碱 基 组 成 


Table 2 Base composition of the mitochondrial AT-rich regions of 13 geographical populations 


of Parnassius glacialis in this study 





居 群 代码 平均 长 度 (bp) 保守 位 点 


变异 位 点 





. ; i . T% C% A% G% (A+T)% 
Population code Average size Conserved site Variable site 

ZJS 492 434 59 48.2 2.4 48.2 1.2. 96.3 
YTS 490 433 62 49.2 2.7 46.7 1.4 95.9 
TTZ 490 409 85 47.3 2.9 48.4 1.4 95.7 
HBY 493 427 70 47.3 2.6 48.7 1.4 95.9 
LJS 493 457 37 47.7 3.0 48.1 1:2 95.7 
TMS 492 445 51 46.5 3.0 49.2 1.2 95.7 
KYS 490 447 48 48.2 2.9 41.6 1.4 95.7 
XLS 492 420 TI 49.0 2.8 41.0 1:2 95.9 
LYS 489 440 53 48.1 217 48.1 1.2 96.1 
NTS 493 434 62 48.1 2.8 47.9 1.2 95.9 
TS 492 401 95 46.5 2.8 49.4 1:2: 95.9 
SNJ 494 415 82 49.6 3.2 45.5 1.6 95.1 
SS 490 411 87 47.3 2.9 47.8 2.0 95.1 





居 群 代码 信息 见 表 1; 下 同 。Population codes are the same as in Table 1. The same below. 
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2.2 冰 清 绢 蝶 地 理 居 群 的 遗传 多 样 性 

325 头 冰 清 绢 蝶 的 线粒体 基因 AT 富 集 区 序列 
比 对 后 显示 ,几乎 所 有 地 理 居 群 内 都 存在 序列 相同 
的 现象 ,但 没有 发 现 完全 共享 单 倍 型 ,但 小 陀 山 单 倍 
型 XLS 6 和 XLS 11 是 小 陀 山 、 黄 柏 源 、 牛 关山、 神农 
架 、 老 君山 天堂 寨 和 天 目 山 所 共享 单 倍 型 ( 表 3)。 
所 有 序列 经 DnaSP 5.10 软件 分 析 , 共 检测 出 239 个 
单 倍 型 (haplotypes) , 单 倍 型 多 样 性 Hd 为 0. 9971, 
TUO LA VEIS AX (nucleotide polymorphism index) 
m H 0. 02948, Theta (per site) Eta f 7j 0.06594, 
就 地 理 居 群 内 部 而 言 , 小 陇 山 单 倍 型 最 多 ,多 达 23 
个 ,其 次 为 牛头 山 (22) 、 神 农 架 (21) 和 黄柏 源 (20)， 
而 云 台 山 (13 ) 单 倍 型 最 少 ;就 单 倍 型 多 样 性 而 言 ， 
老 君 山 (0.989) .小 陇 山 (0.985 ) 和 天 目 山 (0. 985) 
的 单 倍 型 多 样 性 值 较 高 ,而 云 台山 的 单 倍 型 多 样 性 
值 最 小 (0. 940); SE ETE BR E STEM SL 党 
(0.04199) Æ gy (0. 04107 ) 和 神农 架 (0. 04039) 
m ARK, MAA L (0. 01505) m 值 则 较 小 ( 表 1) 。 
这 些 结果 表明 ,不 同 地 理 居 和 群 之 间 已 经 呈现 一 定 程 
度 的 遗传 分 化 ,并 且 小 陇 山 .黄柏 源 、 牛头 山 和 神农 
架 居 群 可 能 为 本 研究 13 个 地 理 居 群 的 祖先 居 群 。 

本 研究 AMOVA 分 析 结 果 显 示 , 冰 清 绢 蝶 各 居 
群 之 间 Fsi 值 大 多 在 0. 15 ~ 0. 25 之 间 , 均 值 为 







































































0.155。 相 比较 而 言 ,泰山 与 天 目 山 居 群 和 神农 架 与 
云 台 山居 群 Fst 值 较 大 ,分 别 为 0.831 和 0.387, 对 
应 的 Nm 值 分 别 为 0.051 和 0.396 ;而 牛头 山 和 小 陇 
山居 群 . 老 君山 和 牛头 山居 群 以 及 神农 架 和 天 堂 寨 
JERE Fst 均 较 小 ,分 别 为 0.050, 0.066 和 0.069, 其 
对 应 的 Nm 值 分 别 为 4.750, 3.538 和 3.373( 表 4) 。 
KLEIT ARR HH , 单 倍 型 XLS 22, LJS 1 和 
XLS 15 EAF KIHI EF (torso), H. XLS 22, 
LJS 1 与 其 他 各 居 群 之 间 联 系 密切 。 但 各 居 群 间 相 
互 交织 在 一 起 ,并 没有 显示 出 明显 的 分 化 格局 ,因此 
这 里 未 显示 梨 支 图 。 
2.3 冰 清 绢 蝶 地 理 居 群 间 的 遗传 距离 

遗传 距离 分 析 表 明 ( 表 5) , 冰 清 绢 蝶 整 个 群体 
的 平均 遗传 距离 为 0. 031 ,遗传 相似 系数 为 0. 969, 
居 群 间 平 均 遗 传 距离 为 0. 033 , 居 群 内 平均 遗传 距 
离 为 0.027。 琅 于 山 和 云 台山 2 个 居 群 的 居 群 内 遗 
传 距离 最 小 (分 别 为 0.017 和 0. 020) , 居 群 内 部 有 
较 多 的 重复 序列 现象 而 存在 相对 较 少 的 单 倍 型 ;而 
神农 架 和 大 别 山 天 和 堂 塞 居 群 内 遗传 距离 较 大 ( 分别 
为 0.038 和 0.033)。 总 体 来 说 , 居 群 内 的 遗传 距离 
普遍 小 于 居 群 间 ,但 老 君 山居 群 较为 特殊 ,其 内 部 个 
体 间 的 遗传 距离 (0. 022 ) 与 其 他 居 群 之 间 的 距离 大 
体 相当 ,其 至 大 于 居 群 间 的 差异 。 





























































































































R3 冰 清 绢 蝶 居 群 内 和 居 群 间 具 有 相同 AT 富 集 区 序列 的 个 体 列表 


Table 3 The sampling individuals of the same AT-rich sequences within and between Parnassius glacialis populations 








居 和 群 代 码 个 体 数 单 倍 型 数目 S 
W 序列 相同 的 个 体 代码 
Population Number of Number of 
Code of the same sequence 
code individuals haplotypes 
ZJS 25 19 ZJS 1/6/7; ZJS 8/10/11/13/; ZJS 3/LYS 19; ZJS 17/TTZ 3 
YTS 25 13 YTS 6/7/9/10/12/15/16/21; YTS 1/5/12/14/22/20 
NS p: bá TTZ 13/14; TTZ 2/6/10/24; TTZ 15/TMS 4/NTS 20/XLS 6/11/NTS 20/LJS 1/SNJ 2/HBY 4; TTZ 25/ 
NTS 23/HBY 25; TTZ 1/4/SNJ 1 
HBY 25 20 HBY 11/12/14; HBY 1/3/6; HBY 15/17; HBY 4/XLS 6/11/NTS 20/LJS 1/SNJ 2/TTZ 15/TMS 4 
Hs m 26 LJS 2/7/17; LJS 3/6/20; LJS 21/XLS 18/NTS 7; LJS 1/SNJ 2/XLS 6/11/NTS 20/HBY 4/TTZ 15/TMS 
i 4; LJS 14/SS 8 
je js T TMS 1/2/3/5/ ; TMS 12/13/20; TMS 6/7/8/9/10; TMS 4/XLS 6/11/NTS 20/LJS 1/SNJ 2/HBY 4/TTZ 
j 15 
TS 25 16 TS 14/11; TS 6/7/8/9/10; TS 1/2/3; TS 20/22/; TS 23/25 
XLS 25 23 XLS 1/4; XLS 6/11/NTS 20/LJS 1/SNJ 2/HBY 4/TTZ 15/TMS 4; XLS 12/SNJ 3/NTS 17 
LYS 25 18 LYS 1/3/5; LYS 8/10/11/12/13/15; LYS 19/NJ 6; LYS 16/LYG 17 
NTS 25 22 NTS 6/7/10; NTS 8/9; NTS 24/SNJ 16; NTS 20/XLS 6/11/LJS 1/SNJ 2/HBY 4/TTZ 15/TMS 4 
KYS 25 15 KYS 1/2/4/5; KYS 3/13/8S 18; KYS 11/12/15; KYS 24/23; KYS 6/7/9/10; KYS 18/YTS 18 
SN] 5s 21 SNJ 6/11/12; SNJ 3/4/7; SNJ 1/TTZ 1/TTZ 4/LYS 6; SNJ 2/XLS 6/11/NTS 20/LJS 1/HBY 4/TTZ 15/ 
à TMS 4 
SS 25 16 SS 1/2/3/4/5; SS 6/1/9/10; SS 13/14/18/19 


居 群 代码 后 的 数字 示 个 体 编号 。The numerals following the population code indicate the individual no. 
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表 4 冰 清 绢 蝶 两 两 居 群 间 的 遗传 分 化 指数 ( Fst) 和 基因 流 (Nm) 
Table 4 The Fst values and the gene flow ( Nm) for pariwise difference of Parnassius glacialis populations in this study 
CAM Fst-inter ZJS YTS TTZ HBY LJS TMS TS XLS LYS NTS KYS SS SNJ 
ZJS 0. 501 1.238 0.080 . 800 .438 .322 1.247 1.869 2.250 1.833 11.462 1.486 0.643 
YTS 0.498 0.168 0.541 .529 .451 .907 0.832 1.059 1.284 1.313 0.964 0.856 0.396 
TTZ 0.514 0.188 0.316 1. 887 .612 .406 1.630 2.955 1.080 1.313 3.128 2.467 3.373 
HBY 0.531 0.238 0.321 0.117 0. 541 .293 2.915 1.816 0.809 0.964 1.221 1.833 1.602 
LJS 0. 538 0.093 0.147 0.290 0.316 .238 0.892 1.451 1.816 3.538 1.221 1.124 0.404 
TMS 0. 524 0.159 0.216 0.151 0.162 0.168 0.051 2.467 2.250 1.688 2.382 2.438 0.814 
TS 0.512 0.167 0.231 0.133 0.079 0.219 0.831 2.528 1.486 1.869 1.943 2.799 1.332 
XLS 0.526 0.118 0.191 0.078 0.121 0.147 0.092 0.090 2.528 4.750 2.382 2.591 1.562 
LYS 0. 504 0.100 0.163 0.188 0.236 0.121 0.100 0.144 0.090 2.599 2.275 2.065 0.633 
NTS 0.526 0.120 0.160 0.160 0.206 0.066 0.129 0.118 0.050 0.089 2.726 2.108 0.775 
KYS 0. 528 0.146 0.206 0.074 0.170 0.170 0.095 0.114 0.095 0.099 0.084 2.762 1.406 
SS 0.522 0.144 0.226 0.092 0.120 0.182 0.093 0.082 0.088 0.108 0.106 0.083 1.486 
SNJ 0. 531 0.280 0.387 0.069 0.135 0.382 0.235 0.158 0.138 0.283 0.244 0.151 0.144 














Fst-inter; 各 居 群 内 部 的 遗传 分 化 指数 The Fst values of the internal population. 对 角 线 下 方 为 遗传 分 化 指数 ( fst) ,对 角 线 上 方 为 基因 流 ( 每 世代 


迁移 数 , Nm) 。The Fst values of the difference between groups are shown below the diagonal, and the Nm (migrants per generation, Nm) values are 


shown above the diagonal. 


表 5 









































in this study based on Kimura-2 parameter model 


冰 清 绢 蝶 13 个 地 理 居 群 的 遗传 距离 (基于 Kimura 双 参 数 模 型 ) 


Table 5 The genetic distance matrices among 13 geographical populations of Parnassius glacialis 








居 群 代码 组 内 组 间 Inter-population 
Population code Intra-population ZJS YTS TTZ HBY LS TMS TS XLS LYS NIS  KYS SS 
ZJS 0. 027 
YTS 0. 020 0.029 
TTZ 0. 033 0.038 0.029 
HBY 0. 023 0.033 0.033 0.031 
LJS 0. 022 0.021 0.020 0.031 0.026 
TMS 0. 023 0.030 0.028 0.032 0.027 0.021 
TS 0. 025 0.032 0.030 0.034 0.025 0.004 0.026 
XLS 0.031 0.034 0.033 0.034 0.030 0.026 0.030 0.030 
LYS 0.017 0.025 0.023 0.030 0.025 0.016 0.022 0.024 0.026 
NTS 0. 029 0.033 0.03 0.037 0.032 0.024 0.028 0.031 0.033 0.026 
KYS 0. 030 0.035 0.033 0.036 0.033 0.027 0.031 0.032 0.033 0.027 0.034 
SS 0.031 0.037 0.036 0.038 0.033 0.030 0.032 0.033 0.036 0.029 0.036 0.037 
SNJ 0. 038 0.045 0.045 0.038 0.034 0.039 0.038 0.036 0.039 0.036 0.043 0.042 0.041 


2.4 冰 清 绢 蝶 地 理 居 群 的 系统 发 生 分 析 
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13 个 地 理 居 群 
分 序列 .AT 富 
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的 12S rRNA 部 


区 RNA“ 和 tRNA 部 分 序列 拼接 
后 的 总 长 度 为 882 ~892 bp , 比 对 后 删除 AT 富 集 区 
中 不 稳定 高 变 区 后 ,经 拼接 得 到 总 长 为 693 bp 的 数 





据 集 , 共 获得 了 113 


他 1 











条 单 倍 型 序列 。 以 白 绢 蝶 等 其 


1 个 绢 蝶 种 类 作为 外 群 
MH561714 - MH561724) ,构建 系统 发 生 树 。 





( GenBank 登录 号 














系统 发 生 分 析 显示 ,尽管 ML 树 的 自 展 值 和 BI 


树 的 后 验 概率 值 不 够 理想 ,但 二 种 方法 构建 的 ML 
树 和 BI 树 拓扑 结构 基本 一 致 。 








明显 


一 








冰 清 绢 蝶 














13 ^J fh 
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rp AGAS RED LS Az LLURUX FH LUCK PR E 
一 起 成 为 一 个 小 支 系 ,而 老 君 山 oem E .泰山 LE LI 
和 云 台山 居 群 大 体 聚 在 一 起 为 一 个 小 支 系 , 且 两 个 
小 支 系 内 部 存在 部 分 混杂 情况 ,如 TTZ 3 与 老 君 山 
居 群 聚 在 一 起 ,ZJS 12 和 ZJS 14 与 云 台 山 和 泰山 居 
群 聚 在 一 起 ,同样 地 ,KYS 6, KYS 8, KYS 21 和 YTS 
2, YTS 3, YTS 4 5E gilt .紫金 山 等 居 群 关系 较为 
密切 (图 1, 2)。 
2.5 冰 清 绢 蝶 地 理 居 群 的 谱系 发 育 年 代 学 分 析 
谱系 年 代 学 结果 表明 ,本 研究 中 所 列 的 绢 蝶 属 
形成 时 间 约 在 中 新 世 早 期 [ 距 今 20. 24 百 万 年 
(Ma), 95% 置信 区 间 : 距 今 24.2 -15.4 Ma] , 冰 清 
绢 蝶 和 白 绢 蝶 为 近 缘 姊妹 种 ,二 者 的 最 近 共 同 祖先 
TERR. 7-25 7. 52Ma 的 中 新 世 晚 期 开始 分 化 。 冰 清 绢 
蝶 的 形成 时 间 约 为 第 四 纪 更 新 世 早 期 ( 距 今 1.53 
Ma, 9596 "Efi EX JR]: IES 2. 54 -0.51 Ma) ,第 工大 
支 系 各 群体 形成 和 群体 间 分 歧 时 间 大 多 在 更 新 世 早 
WOES 1.53 -1.14 Ma) , 28 ILA SC SJÉ ROT ATE 
的 时 间 主 要 在 更 新 世 中 晚期 (0. 81 -0.28 Ma) ,两 
大 支 系 扩张 交流 的 时 间 主 要 在 更 新 世 中 期 ( 距 今 
1.14 -0.81 Ma) 。 冰 清 绢 蝶 的 形成 时 间 、 两 大 文系 
的 分 化 时 间 与 扩张 交流 时 间 都 在 第 四 纪 冰 期 中 的 冰 
期 间 冰 期 的 轮回 事件 之 内 (图 2) 。 
2.6 冰 清 绢 蝶 扩 散 路 径 分 析 
依据 现今 中 国 冰 清 绢 蝶 的 地 理 分 布 ,结合 本 研 
究 冰 清 绢 蝶 各 居 群 的 谱系 演化 分歧 时 间 及 其 遗传 
分 化 状况 的 结果 , 冰 清 绢 蝶 祖 先 可 能 于 我 国 原 西南 
地 区 ( 现 青藏 高 原 东北 缘 及 横断 山脉 一 带 )( 汪 学 从 
等 , 2013; 张 亚 玲 等 , 2015; 王 运 良 , 2017) ,在 距 今 
约 1.53 Ma 的 渐 新 世 时 期 开始 形成 和 分 化 ,之 后 沿 
着 秦岭 山系 进行 扩散 。 随 着 第 四 纪 的 冰期 - 间 冰 期 
轮回 事件 , 冰 清 绢 蝶 从 秦岭 山系 向 东北 、 东 南 形成 2 
大 文系 迁移 到 大 别 山 系 和 伏 牛 山系 ,并 沿 着 这 些 山 
系 继续 向 南 或 向 东 迁 移 ,向 南 甚 至 到 达到 宁 镇 山脉 
(天 目 山 、 紫 金山 ) ,向 东 到 达 泰 山 山脉 (泰山 、 昆 寄 
山 云 台山 等 处 ) (图 3) 。 


































































































3 讨论 





昆虫 线粒体 基因 组 AT 富 集 区 属于 不 编码 基 
,在 进化 过 程 中 所 承受 的 选择 压力 较 小 ,同时 ,和 
线粒体 基因 组 其 他 部 分 相 比 ,其 进化 速率 较 快 ,种 间 
甚至 种 内 间 变 异 明显 , 因 此 ,它们 在 亲缘 关系 较 近 的 
种 间 系 统 分 类 种群 分 化 和 种 内 遗传 多 样 性 等 研究 





























中 作为 有 效 的 分 子 标记 而 被 广泛 运用 ( Wolstenholme , 
1992; Hao et al., 2013; 王 维 等 , 2013; Kim et al., 
2014; 何 海燕 等 , 2016)。 本 研究 表明 , 冰 清 绢 蝶 的 
线粒体 AT 富 集 区 序列 具有 绝 大 多 数 鳞 怒 目的 典型 
村 征 ,如 ATAGA (T) 引导 的 Ploy T 结构 、(AT)n 或 
35 (TA) n 元 件 以 及 3' 端 的 Poly ( A) 结构 等 (Hao et 
al., 2013; 王 维 等 , 2013; Chen et al., 2014; Wang 
et al., 2015) 。 虽 然 其 序列 长 度 相 对 较 短 ,但 在 居 群 
间 和 居 群 内 仍 存 在 较为 丰富 的 多 态 性 ,是 探讨 冰 清 
绢 蝶 群 体 遗 传 分 化 和 谱系 发 生 较 为 理想 的 分 子 标记 
之 二 

本 研究 显示 冰 清 绢 蝶 的 单 倍 型 多 样 性 Ha 为 
0. 9971 , 核 背 酸 多 态 性 指数 m 值 为 0.02948( 表 1)。 
依据 Grant 和 Bowen(1998) 等 的 观点 , 即 Hd z0.5, 
m=0. 005 时 ,表示 种 群 进化 历史 较为 悠久 ,种群 较 
为 稳定 。 因 而 , 冰 清 绢 蝶 在 我 国 可 能 经 过 了 较为 漫 
长 的 进化 历史 , 才 得 以 形成 当今 的 地 理 分 布 格局 。 
另外 , 本 人 研究 估算 的 冰 清 绢 蝶 的 平均 Fst 值 为 
0. 155 ,依据 Wright (1990). 提出 的 观点 , BH Fst 为 
0. 15 ~0.25 时 ,说 明 群 体 间 遗 传 分 化 相对 较 大 , 则 
中 国 冰 清 绢 蝶 有 一 定 的 遗传 分 化 ,其 遗传 分 化 主要 
来 自居 群 内 部 。 本 研究 13 个 居 群 的 Nm 平均 值 为 
1.50, 绝 大 部 分 两 两 居 群 间 的 基因 流 Nm 值 都 >1， 
表明 居 群 间 存 在 一 定 的 基因 交流 。 值 得 注意 的 是 小 
陇 山 和 牛头 山居 群 间 以 及 老 君 山 和 牛头 山居 群 间 的 
Nm 值 分 别 为 4.750 和 3. 538 ,表明 这 些 居 群 间 存在 
较为 明显 的 基因 交流 ,应 与 秦岭 山系 和 伏 牛 山系 间 
近 距 离 扩 散 有 关 ; 同 样 ,神农 架 和 天 堂 寨 居 和 群 间 也 存 
在 类 似 的 现象 (Nm 值 为 3.373) ,可 能 与 秦岭 山系 和 
大 别 山 系 相互 交流 和 扩散 有 关 。 

冰 清 绢 蝶 居 群 间 的 平均 遗传 距离 (0. 033) 略 大 
于 居 群 内 平均 遗传 距离 (0.031)。 但 老 君 山居 群 内 

异 以 及 该 居 群 与 其 他 居 群 间 的 遗传 距离 没有 明显 

差异 ,特别 是 与 小 陇 山 .牛头 山 嵩山 .泰山 等 居 群 间 
的 差异 其 至 小 于 其 居 群 内 部 差异 ,表明 其 遗传 分 化 
多 在 居 群 内 部 ( 表 4) 。 同 时 ,本 研究 的 梨 枝 图 (未 显 
示 ) 也 表明 各 居 群 间 相 互 交织 在 一 起 ,没有 明显 的 
分 化 格局 。 最 后 ,从 系统 发 生 树 (图 1, 2) 来 看 , 尽 
管 老 君 山居 群 的 LJS 9, LJS 17 和 LJS 25 单 倍 型 与 
第 工 文系 相互 聚 在 一 起 , 而 黄柏 源 居 群 的 HBY 4, 
HBY 14 和 牛头 山居 群 的 NTS 10 和 NTS 15 以 及 神 
农 架 的 SNJ 18 与 第 开支 系 聚 集 在 一 起 ,但 总 体 来 
看 , 冰 清 绢 蝶 大 体 上 仍 分 化 为 2 大 支 系 。 

上 述 结果 提示 ,在 第 四 纪 冰 期 - 间 冰 期 的 轮回 事 












































































































































































































































王 运 良 等 : 


483 


基于 线粒体 基因 组 控 表 





站 


区 序列 的 中 国 冰 清 绢 蝶 遗 传 分 化 和 谱系 生物 地 理 研究 





| 




















illt Parnassius c 
It Parnassius a: 








LYS 6 
LYS 
ZJS6 
ZJS 10-ZJS9 
ZJS2 
LYSS 
XLS 12 x 
XLS 16* 
SS 18 


I 


SS 14 
KYS 





s510 —— 

TMS 19-TTZ 13-TTZ 9 
TMS 18 

TMS 14 


SS 17 
SS 18 
Kilometers 





HBY 25-NTS 23 
HBY 13 
SUE 
SNJ 15 
SNI2 


HbY 23 
LJS22 x 
LJS 23 * 
HBY 11-HBY 5 
SNJ 25 
XLS 21 
SNJ 3-SNJ 4-NTS 17 


XLS 12 

XLS 16-XLS 18-NTS 7 

XLS 24-HBY 11-HBY 14 
TS 25 


XLS 19 
XLS 24 
SNJ 16 
XLS3 
NTS 1 
NIS 2-NTS 5 
BY 8-HBY9 
XLS 24 
LJS 3-LJS 2 
LJS 2-LJS 17-LJS «x 
SNJ 14 


SNJ 17 
HBY 7-NTS 1 
HBY 17-HBY 20 
SNJ 20 
XLS 18 
LJS21 > 
HBY 3 
HBY 2 
HBY 2 
HBY 18-HBY 24 
LIS 25 > 
SNJ 2-SNJ 8 
XLS 20-NTS 20 
HBY 4-HBY 5-HBY 7 
XLS 11-X1S 23 
XIS 11 








Ht Parnassius stubbendorfii ( MH561716) 
BSR Parnassius orleans (MH561717) 
度 绢 蝶 Parnassius andrejii (MH561714) 
mo (MH561719) 
HE Parnassius imperator (MH561720) 
红 : € Parnassius bremeri (MH561723) 
mE AE Parnassius apollo (MH561722) 
依 帕 绢 蝶 Parnassius epaphus (MH561721) 
SERI Parnassius apolloinus (MH561724) 
:ephalus (MH561715) 
cco (MH561714) A 


0.04 








图 1 基于 Data 1 构建 中 国 冰 清 绢 蝶 13 个 地 理 居 和 群 ML 系统 发 生 树 








Fig. 1 The ML tree of the 13 geographic populations of Parnassius glacialis in China based on Data 1 
A: JEW , 星 号 示 集 群 中 交叉 的 单 倍 型 Rectangular tree, stars indicating the hylotypes affiliated to different lineages; B; 放射 树 Radical tree; C; I 



































和 开 两 大 集群 在 地 图 上 位 置 ， 
山 \ 昆 帘 山 和 云 台 山居 群 Two li 
I includes the populations TTZ, 
部 分 序列 构成 的 序列 集合 。 地 
号 个 体 的 序列 代码 ,其 余 类 推 。 


sequence and tRNA" gene partia 








is the sequence code of the No. 1 


for Fig. 2. 


中 工 主要 为 黄柏 源 、 小 陇 山 .牛头 山 和 神农 架 居 群 ， 开 主要 为 天 堂 寨 、 紫 金山 RR EL LU .嵩山 gA R 
neages ( [ and [[ ) on the map. The lineage I includes the populations HBY, XLS, NTS and SNJ, while the lineage 
ZJS, LYS, TMS, SS, LJS, TS, KYS and YTS. Data 1 为 128 rRNA 部 分 序列 AT 富 集 区 序列 RNAMe 和 tRNAL 
图 文件 信息 来 自 国家 基础 地 理 信 息 中 心 (http: // ngec. sbsm. gov. cen/) o XLS 11 为 冰 清 绢 蝶 天 水 小 陇 山居 群 11 
Data 1 is the dataset of combined 12S rRNA gene partial sequence, AT-rich region sequence, tRNAM® gene partial 
sequence. The map information is from the National Geomatics Center of China ( http: // ngcc. sbsm. gov. cn/). XLS 11 
1 individual of the Xiaolongshan population of P. glacialis, and the others may be deduced by analogy. 图 2 [i] The same 
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图 2 基于 Data 1 构建 中 国 冰 清 绢 蝶 13 个 地 理 居 群 BI 系统 发 生 树 和 谱系 发 育 年 代 图 
Fig. 2 The BI tree and phylochronological chart of 13 geographic populations of Parnassius glacialis in China based on Data 1 
A [pos OL 2 E E 8], CD, @ 和 图 箭头 所 示 为 分 子 钟 校正 点 位 置 。The nodal values show the divergence time in Fig. A. (D, © and © 


represent the calibration points of molecular clock. 














(EF, RKD TRECE, EEKMA —— 起源 和 分 化 的 中 心地 区 。 另 外 ,本 研究 显示 ,河南 老 
架 等 地 的 冰 清 绢 灵 居 和 群 基因 交流 明显 , 且 它 们 在 聚 ”君山 居 群 的 部 分 单 倍 型 在 系统 树 的 2 个 支 系 中 都 有 
类 中 并 没有 很 好 区 分 开 , 因 而 ,这 一 地 区 是 冰 清 绢 蝶 ”分 布 ,因而 老 君 山居 群 可 能 是 我 国 冰 清 绢 蝶 发 生 频 
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图 3 中 国 冰 清 绢 蝶 扩 散 路 径 示 意图 
Fig. 3 Schematic diagram of the dispersal route of Parnassius glacialis in China 


地 图 信息 来 自 国家 基础 地 理 信息 中 心 (http: // ngec. sbsm. gov. cn) The map information is from the National Geomatics Center of China (http: // 



































ngcc. sbsm. gov. cn). 


繁 交 流 的 交汇 中 心 之 一 。 再 者 ,虽然 在 后 续 的 分 化 
过 程 中 ,大 别 山 天 党 寨 居 群 和 琅 志 山 、 天 目 山 、 紫 金 
山居 群 构成 一 个 小 的 支 系 ,而 老 君 山居 群 与 党 山 、 泰 
山 \ 云 台山 、 昆 沉 山 等 居 群 构成 男 一 个 小 支 系 。 但 上 
述 2 个 小 支 系 的 居 群 间 均 存在 部 分 相同 的 单 倍 型 ， 
尤其 是 位 于 江苏 连云港 的 云 台 山居 群 ,尽管 其 单 倍 
型 较 少 ,但 很 明显 的 与 两 支 系 都 有 交叉 ,在 老 君 山 等 
居 群 构成 的 小 支 系 中 , 云 台 山居 和 群 与 昆 咀 山居 群 聚 
在 一 起 ,而 在 大 别 山 天 党 寨 等 居 群 中 构成 的 男 一 文 
小 文系 中 , 云 台 山居 群 与 紫金 山 ` 天 目 山 居 群 亲缘 关 
系 较 近 ,因而 , 云 台 山居 群 应 与 紫金 山居 和 群 和 昆 般 山 
等 居 群 也 发 生 了 某 种 程度 的 基因 交流 。 

第 四 纪 以 全 球 气候 显著 变 冷 为 其 重要 特点 , 且 
形成 多 次 冰期 - 间 冰 期 的 轮回 ,这 种 气候 剧烈 波动 ， 
引起 了 生物 大 规模 的 迁移 ,甚至 是 生存 灾难 (Andel 
and Tzedakis, 1996; 陈 冬 梅 等 , 2011) 。 前 期 的 分 子 
钟 研究 标明 , 绢 蝶 属 可 能 起 源 于 古 近 纪 晚 期 或 新 近 
纪 早 期 ( 距 今 35 - 17 Ma) 的 我 国 西南 地 区 ,现今 的 
横断 山脉 一 带 ( 汪 学 俭 等 , 2013; 张 亚 玲 等 , 2015; 
EZR, 2017) ,中 新 世 ( 距 今 23 -5 Ma) 以 来 ,我 国 
的 喜马拉雅 造山 运动 加 剧 ,青藏 高 原 快 速 抬 升 ,天 
山 ,秦岭 等 山系 也 继续 抬升 , 绢 蝶 并 产生 了 一 些 适 应 































































































新 世 后 期 , 冰 清 绢 蝶 和 其 姊妹 种 白 绢 蝶 的 共同 祖先 
开始 发 生 分 化 , 距 今 约 1.7 Ma 的 渐 新 世 早 期 , 冰 清 
绢 蝶 这 一 物种 最 终 形成 并 开始 分 化 。 

随后 而 来 的 第 四 纪 ( 距 今 2.6 Ma - ) 冰 期 - 间 冰 
期 轮回 事件 .寄主 植 物 的 扩散 及 其 他 地 球 环境 事件 
( 孙 航 ，2002 ; Zhang et al., 2016) ,9| 起 冰 清 绢 蝶 不 
断 地 扩散 、 隔 离 以 及 种 群 间 的 相互 混杂 现象 。 特 别 
是 , 当 冰 期 来 临时 ,其 寄主 植物 ( 延 胡 索 Corydalis 
yanhusuo 或 黄 曹 Corydalis pallida ) 向 南方 或 向 低 海 
拔 的 地 方 退缩 ; 当 间 冰期 期 间 ,海平 面 上 升 ,温度 升 
高 ,植物 再 向 北方 或 高 海拔 地 方 扩散 ,而 绢 蝶 和 其 他 
生物 一 样 , 随 着 其 寄主 植物 一 起 来 回迁 移 。 结 合 本 
研究 冰 清 绢 蝶 各 居 群 的 聚 类 情况 分歧 时 间 及 其 分 
布 来 看 , 冰 清 绢 蝶 祖 先 可 能 先 从 横断 山脉 附近 经 青 
藏 高 原 东北 缘 扩 散 到 秦岭 一 带 , 之 后 沿 东 北 和 东南 
向 二 条 主要 路 径 由 秦岭 地 区 逐渐 向 其 他 低 海 拔 地 区 
扩散 (其 中 LB 29 1.8 Ma, 随 着 第 四 纪 冰 期 的 不 断 
发 展 , 因 海 平面 下 降 和 陆桥 的 形成 等 事件 而 由 中 国 
大 陆 沿 海区 域 扩散 到 日 本 地 区 ,并 随后 由 于 日 本 海 
的 渐进 扩张 而 导致 的 地 理 隔离 逐渐 分 化 为 一 个 新 的 
亚 种 ( Yagi et al., 2001)) 。 

随 着 第 四 纪 冰 期 的 持续 发 展 ,尤其 是 在 西夏 邦 



































高 海拔 环境 的 特性 ,如 翅 被 毛 .耐寒 、 耐 低 氧 等 特征 。 
之 后 , 随 着 栖息 环境 的 快速 变化 和 种 间 苋 争 的 加 剧 ， 
绢 蝶 属 发 生 快 速 适应 性 分 化 ;在 距 今 约 7.5 Ma 的 中 








马 冰 期 (Shishapangma Glaciation) ( HES 1. 17 -0.8 
Ma) ,温度 持续 下 降 , 寄 主 植物 可 能 在 这 一 时 期 得 以 
向 低 纬度 或 低 海拔 地 区 扩张 , 冰 清 绢 蝶 也 得 以 从 秦 
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岭 山 系 向 大 别 山系 和 伏 牛 山系 扩散 ,这 一 时 间 段 恰 
好 与 冰 清 绢 蝶 第 1 大 文系 和 第 2 文系 扩张 交流 的 时 
间 ( 距 今 1. 14 -0.81 Ma) 相 吻 合 而 得 以 相互 验证 。 
随后 , 随 着 望 坤 冰期 (Wangkun Glaciation ) ( JE S 
0.76 -0.5 Ma) 地 球 温度 的 急剧 降低 ,冰川 面积 扩 
大 ,高 纬度 或 高 山高 原 的 生物 继续 向 低 纬 度 、 低 海拔 
地 区 或 向 南 迁 移 , 望 坤 冰期 后 , 随 之 而 来 的 是 第 四 纪 
冰期 中 时 间 最 长 的 一 次 间 冰 期 ,持续 时 间 近 100 
Ka, 这 期 间 , 由 于 温度 升 高 ,海平 面 上 升 , 冰 清 绢 蝶 
开始 向 高 纬度 或 高 海拔 地 区 迁移 ,也 有 一 部 分 遗留 
在 低 海 拔 的 高 山地 区 (如 天 目 山 、 云 台山 \ 昆 常山 、 
泰山 等 地 ) ,同时 ,这 次 间 冰 期 可 能 促使 了 伏 牛 山系 
和 大 别 山系 的 冰 清 绢 蝶 进 行 了 基因 交流 。 随 后 义 伴 
随 着 多 次 冰期 - 间 冰 期 轮回 和 东南 季风 事件 (Zhang 
et al., 2016) , 宁 镇 山脉 的 种 群 也 和 泰山 山脉 的 种 篆 
也 发 生 了 进一步 交流 。 

综合 多 方面 证 据 分 析 , 冰 清 绢 蝶 不 同居 群 间 的 
扩散 路 径 与 山 肪 分布. 走向 以 及 地 理 隔离 等 存在 明 
的 相关 性 ,并 受到 第 四 纪 至 今 的 气候 变化 .生物 环 
境 等 因素 的 影响 , 才 最 终 形成 了 今天 的 冰 清 绢 蝶 分 
布 格局 。 
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